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Abstract 

We derive a simple set of non-linear, (1 -|- l)-dimensional partial differential equations 
that describe the dynamical evolution of black strings and branes to leading order in 
the expansion in the inverse of the number of dimensions D. These equations are easily 
solved numerically. Their solution shows that thin enough black strings are unstable to 
developing inhomogeneities along their length, and at late times they asymptote to stable 
non-uniform black strings. This proves an earlier conjecture about the endpoint of the 
instability of black strings in a large enough number of dimensions. If the initial black 
string is very thin, the final conhguration is highly non-uniform and resembles a periodic 
array of localized black holes joined by short necks. We also present the equations that 
describe the non-linear dynamics of Anti-deSitter black branes at large D. 



I. The instability of black strings and black branes discovered in [1] is a phenomenon with 
wide implications for the physics of higher-dimensional black holes and their applications 
to string theory and gauge/gravity duality [2]. Black strings behave similarly to tubular 
soap films and are prone, when thin enough, to rippling, i.e., developing non-uniformities 
along their direction. This instability is well established in perturbation theory, but its 
growth beyond the linearized approximation and its endpoint at asymptotically late times 
are notoriously hard problems. The numerical evolution of a perturbed five-dimensional 
black string in [3] is a landmark result: it gives strong evidence that the classical evolution 
does not stop at any stable configuration but proceeds in a self-similar cascade to smaller 
scales. 

Important as the result of [3] is, it still leaves open many questions about the fate 
of unstable black strings and black branes. For instance, it would be very convenient to 
have a better, possibly analytic, understanding of the late-time evolution. Moreover, [3] 
presents a single calculation of a single system. Do all unstable black strings behave in the 
same manner? This is indeed unlikely: there has long been evidence that, by modifying 
the parameters of the system, the endpoint may be different. In the simplest instance, the 
parameter is the number of spacetime dimensions where the black string lives. Ref. [4] 
found a critical dimension, D = D* ~ 13.5, above which weakly non-uniform static black 
strings (NUBS) have larger horizon area, for fixed string length and mass, than the uniform 
solutions. It is then possible that the classical evolution of the system in D > D* ends 
at a stable NUBS, as proposed (independently of D) in [5]. However, given the almost 
prohibitive cost of the numerical simulations (100 000 CPU hours for [3]), the investigation 
of this possibility and the systematic study of the evolution of related systems has not 
been undertaken to this date. This is strong motivation to search for simpler methods 
capable of capturing at least the main qualitative features of the phenomenon. 

In this article we present an approach, based on an expansion in 1/D, that allows to 
address some of these issues. In particular, we show that for large enough D the end¬ 
point of the instability is generically a stable NUBS. The simplification of the problem 
is dramatic: accurate numerical evolutions can be obtained in seconds (or less) in a con¬ 
ventional computer running a one-line NDSolve of Mathematica. This encourages further 
investigation of this and similar problems by these means. 

II. The large-D approach to black hole physics, initiated in [6, 7, 8], has been recently 
developed to deal with fully non-linear deformations of horizons [9, 10, 11, 12]. While the 
formalism of [10] incorporates time evolution, it is unclear to us whether it allows, as it 
is, to study the system at hand, which involves horizon length scales 1/^/D. Therefore, 
we solve the problem ab initio in a formulation adapted to the dynamics of black branes. 

Using ingoing Eddington-Finkelstein coordinates, the metric of a uniform black p- 
brane in D = n + p + 3 dimensions, boosted along its world volume with velocity vector 
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Ufj,, r]f^'^Ufj,Uu = -1, = 0,... ,p, is 

ds^ = u^Uu'j da^da" - 2u^da^dr + r'^dVLn+i • (1) 

The worldvolume directions are 


cj^ = (t,o-“), a = l,...,p. 


( 2 ) 


We take p to be finite, so -D —>• oo is n —)• oo, and we will study deformations of the black 
brane that depend on a^. 

Since the wavelength of unstable fluctuations is known to be ~ \fn [6, 7], we rescale 
(t“ —)• cj'^!y/n. Furthermore, we consider small velocities, = (—1 + O (n“^) ,Valy/n) , 
so, denoting m = Tq , the metric (1) becomes 




- (l - ^) df + 2(^dt- ^da^-) dr - 

+- (dab + da^-da’’ + r'^dQn+i 

n \ nr^ / 


2mva 
n r”' 


da°'dt 


(3) 


We seek solutions that can be regarded as having m and Va not constant but varying 
with the coordinates In order to find them, we take a metric ansatz of the form 


ds^ = —Adt^ — 2{utdt + Uada°‘)dr — 2Cada°‘dt + Gabda^da^ + r‘^dQn+i, (4) 


with 


■^ = E h • 




a^, r) 


k>0 


k>0 
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u'a\cr^,r) 


Ca = Y. 




fc >0 
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Gab=-{ dab + V 
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(5) 


( 6 ) 


The different scalings with n conform to (3). We have chosen a Bondi-type gauge where 
grr = 0 and r is the area-radius of 5'"'"’“^ to all orders in 1/n. This leaves a gauge freedom 
in the choice of Ua, which is a shift vector on surfaces at constant r and is only restricted 
by boundary conditions. We have partially gauge-fixed it to be independent of r. 

We introduce the radial coordinate 


R= ( - 

To 


(7) 


such that when D ^ oo keeping R hnite, we focus on the near-horizon region. The 
boundary conditions at large R (for near-horizon decoupled geometries [13]) are 
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T = l + 0(R-i), Ca = 0{R-^), Gab = -{dab + 0{n-\R-^)) 


( 8 ) 


We also require regularity at the horizon, which is where g^i^u^u’^ = 0. Henceforth we fix 
ro = 1. 
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It is now straightforward to solve the Einstein equations perturbatively in 1/n. To 
leading order, the solution is 

^( 0 ) ^ ^(0) _ MO) _ Pa{(y^)Pb{(y^) 


= 


^( 0 ) ^ Pa{(y^)Pb{<y^) 
m{a^^)R 




= constant. 


We have gauge-fixed the cr^-dependence of Ua in a manner consistent with boundary 
conditions. Solving similarly at the next order, the equations Ru = 0 and Rta = 0 are 
r-independent constraints that require, respectively. 


dtm — dbd^m = —dbp^ , 


and ^ 

dtPa - dbd'^pa = dam - db (12) 

\ m J 

{a,b indices are raised with the flat metric 5“^). These equations are one of our main 
results. They are the effective equations for the collective variables Pa{cr^) that 

describe non-linear fluctuations of the black brane. These variables give the energy and 
momentum densities of the black brane on sections at constant M at a given time. The 
horizon of the solution is at R = so the function m{a^) is also interpreted as 

the area-density, i.e.,the area of the (up to a factor of the unit-sphere area) at the 

horizon at a fixed value of t and cj“, to leading order in 1/n. For solutions that are spatially 
periodic in the quantities 


M = J da^ ... da^ , Pa = J da^ ... da^ pa{t,a^) 


are conserved in time. Up to normalization factors, they are the total mass and momentum 
of the black brane. Up to similar factors and to leading order, M is also the total horizon 
area, which therefore does not vary as the system evolves. If we define the surface gravity 
as the non-affinity of the horizon generator dt, it also remains constant. 

The equations are invariant under Galilean boosts along cj“ with constant velocity Va, 

M^M-Vat, Pa^Pa + mVa, (14) 

which allows to fix the rest frame of the black brane, in which Pa = 0, and set = 0. 

This effective formulation of black brane dynamics passes two important checks: (i) 
For small, linearized perturbations around the uniform black brane, with momentum k 
aligned with = z, 

m{t, z) = l + Sm , p,(t, z) = 6pa , (15) 

the solution frequencies are 

uj± = i{±k — k ^), u = —ik^ , (16) 
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which reproduce the frequencies of the sound and shear modes of the black brane to 
leading order at large D [7]. In particular, the frequency a;_|_ for 0 < A; < 1 corresponds to 
the Gregory-Laflamme unstable mode, (ii) If we consider static, shear-free deformations, 
m = m{z), Pi = p{z), PajLi = 0, the resulting equation (with p = m' in the rest frame) 


m +m — 


{m'y 


= constant. 


(17) 


m 

is equivalent to the equation for static black strings derived in [9, 12] (using different 
gauges), with the variables being related by m{z) = exp(2P(z)). 

Observe that the sound deformations m, pi along a direction cr^, are not affected 
by the shear Pa^i- In the following we set Pa^i = 0 and consider black strings with 
pi{t,(7^ = z) = p{t,z). The direction is compactified, z G [—L/2,L/2]. We parametrize 
the periodicity L in terms of a wavenumber kL as 

tL = f . ( 18 ) 


Since we are fixing the string thickness rg = 1, the uniform strings are characterized by 
the value of ki- Smaller values of k^ correspond to thinner black strings. 


III. We solve numerically the black string equations 

dtm{t, z) - z) = -dzp{t, z ), (19) 


and 2 

dtp{t, z) - dlp{t, z) = dz z) - ^ . (20) 

The Mathematica function NDSolve handles them without difficulty. We fix a value of 
and introduce a small perturbation of the static uniform black string, m(0, z) = l+6mo{z), 
p{0,z) = 6po{z), such that the momentum P vanishes. We find that for kL > 1, the 
perturbation quickly dissipates and the black string becomes uniform, in agreement with 
the absence in (16) of unstable linear modes with wavenumber A: > 1. 

For thinner black strings, with A:^ < 1, after brief initial transients the deformation 
grows (at approximately the exponential rate of the linearized solution). Eventually, the 
system settles down at a stable configuration that approximates very well a solution of the 
static equation (17). M and P in (13) remain constant to very good accuracy throughout 
the evolution. 

In figures 1 and 2 we show two sample simulations. We plot m(t,z), which gives the 
area of the at a given z. This is different than the area-radius of these spheres, which 
is 


TZ{t, z) = m{z) »i+i ~ 1 -|- 


lnm(t, z) 


n 


( 21 ) 


For large deformations. In m{t,z) can reveal structure that in m{t,z) is exponentially 
suppressed (insets in fig. 2). 
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Figure 1: Dynamical evolution of a perturbed string with = 0.98. The horizontal axis is the 
^n+i_area function m{t,z). For this simulation, the hnal state was reached before t = 500. The 
time it takes depends on the size and shape of the initial perturbation, but the final configuration 
is independent of them. 


Fig. 1 is the evolution of a not-too-thin black string with = 0.98. Since the per¬ 
turbative unstable rate is small, the evolution is slow. The final profile is approximately 
sinusoidal, as expected for a small deformation. 



Figure 2: Dynamical evolution of a perturbed string with = 0.55. For the last two plots we 
display him{t, z) in insets, which show that the blobs extend to almost fill the compact direction. 
For this simulation, the final state was reached around f = 40. 


Fig. 2 shows the evolution for a thinner black string with kL = 0.55, which evolves 
much faster and develops a large blob in its hnal state. As shown in [9, 12], the prohle of 
large static blobs is very approximately gaussian, 

m{z) ~ (22) 

(the amplitude is hxed to have the same M as & uniform black string of length L). The 
radius for (22) decays only like — 2^/2 away from the center, as seen in the insets in 
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fig. 2, and the blob fills almost all the compact direction. Indeed, (22) gives an excellent 
approximation to the profile of Inm(z) for the final solution at all z (better than 1% for 
ki < 0.6) except in a very short ‘neck’ of length (A2;)neck — {^/L) Ini near \z\ = L/2 [9]. 
As shown in [12], these blobs approximate very well the shape of a spherical black hole 
within a region of width ~ ^I'/n around its equator, as well as its total area. Bear in 
mind, however, that while the height of the blobs is TZ = 0(1), their proper length along 
z is L/y/n, i.e., they are much thinner than a sphere. 

When 1/2 < /c£, < 1 the NUBS at the endpoint of the evolution is unique for each 
value of ki, independently of the shape and amplitude of the initial perturbation. When 
0 < /ci < 1/2 there are (at least) two different unstable modes that can be excited, with 
wavenumbers /cl and 2kL (or a higher multiple), and two possible final static NUBS. Since 
their mass and area are the same to leading order in 1/n, the evolution depends on the 
relative growth rates of the unstable modes and on the shape of the initial perturbation. 
When the solution develops two or more blobs, it may be difficult to determine if these will 
remain in the asymptotic final state, or are instead only part of a very long-lived transient 
phase: since m is exponentially small inbetween the blobs, the interaction among them 
may be lost in numerical error. These situations, however, may fall outside the range of 
applicability of this formulation. 

The large-T> approximation is valid when jlnmj, Inmj, Inpj <C n. In the 
solutions we consider, the most stringent condition is jlnmj <C n, which at the neck 
becomes, L <C 2y/^, i.e., nk\ S> 7r^/2. Thus our results are quantitatively accurate 
in a given dimension only for sufficiently large k^, or conversely, for a given k^ only in 
sufficiently large re. However, it may be reasonable to expect that the model still captures 
correctly the qualitative physics for re above the critical dimension re* = Z?* — 4 ~ 9.5. 
Observe that near the limit of validity, at L ~ y/n, the blob has proper size 0(1) in all 
directions [9]. This suggests that this final state can be regarded as an array of roughly 
spherical black holes joined by thin, short necks. Additionally, when the neck becomes thin 
in Planck units, the black string may break up into separate black holes due to quantum 
gravity effects. 

Within this range of validity, we conclude that our results are very strong evidence 
that the endpoint of the black string instability at large enough D is generically a stable 
NUBS. 

IV. There are several possible extensions of our study. With little extra effort we can 
obtain the equations that describe the leading large-D non-linear dynamics of Anti-deSitter 
black branes. Either by directly solving the equations as in the previous case, or by 
applying the AdS/Ricci-flat correspondence of jl4], we find 

dtm — dbd^m = —dtp^ , (23) 
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and 



dtPa - dbd'^pa = -dam - db 


(24) 


It is easy to prove that these equations do not have any non-uniform static solutions. Small 


dynamical perturbations of the uniform configuration, (15), give as solutions the sound 
and shear quasinormal frequencies of the AdS black brane computed in [15], namely, 


ui± = ±k — ik^ , oj = —ik^ . 


(25) 


Other quasinormal modes of the black brane have frequencies 0{D) and do not appear 
in the near-horizon decoupling limit. Allowing variation in a number p ~ Z) of horizon 
directions may involve qualitative changes. 

The inclusion of 1/D corrections to the equations (11), (12) is potentially very inter¬ 
esting. It has been shown in [12] that 1/D corrections allow to accurately identify the 
critical dimension D* below which static NUBS of a given area have higher mass than 
uniform black strings, so they cannot be the endpoints of the instability [16, 17]. This 
suggests that the 1/D expansion may be able to reproduce, when D < D^, an evolution 
qualitatively similar to the one observed in [3]. We hope to report on this in the future. 
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